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Abstract 

Einstein's equations admit solutions corresponding to photon rockets. 
In these a massive particle recoils because of the anisotropic emission of 
(•^ ' photons. In this paper we ask whether rocket motion can be powered 

("^ , only by the emission of gravitational waves. We use the double series 

^S) ' approximation method and show that this is possible. A loss of mass 

c7^ ' and gain in momentum arise in the second approximation because of the 

C^^ . emission of quadrupole and octupole waves. 

Q-i! 1 Introduction 

5— I . 

f^j\, Recently there have appeared a number of investigations into exact solutions of 

Einstein's equations corresponding to photon rockets [l]-[7]. In these a massive 

particle recoils because of the anisotropic emission of photons. One can ask 

^^ ' whether a similar motion is possible powered by the emission of gravitational 

waves rather than photons. An answer to this question is probably implicit 
in work many years ago by Newman and collaborators [8]- [12] who studied 
the motion of a gravitating system subject to no external forces but taking 
into account its internal structure. Another related work is that of Hogan and 
Robinson [13] on the deviation of a particle from geodesic motion under the 
influence of its own gravitational radiation. However these workers were not 
concerned explicitly with rocket motion, and it is not clear to us how their 
results could be applied to the problem. 

In this paper we study the gravitational wave rocket, using an approx- 
imation method. We start with the photon rocket metric 



ds^ 


= (1- 


- fr^ sin^ 9 - 2/rcos6' - 2Mr-^)du^ + 2dudr 
-r'^idO'^ + sin^ ed4)'^) - 2fr'^ sin OdudO 



(1) 



with 

< r, <9 <TT, < (J) <2n, -oo < u < oo. 

The coordinates will be numbered 

X — r, X = 6, X — (p, X — u. 

M and / are arbitrary functions of the retarded time u, to be interpreted as the 
mass of a particle at the origin and its acceleration along the negative direction 
of the polar axis. If M = the spacetime is flat and the coordinate system has 
an origin with acceleration /. Metric (1) satisfies Einstein's equations 

R''' - (l/2)i?g*'= = -SnT'' = -STTfiPl'' (2) 

where /i is a scalar function and P is a null vector. The right-hand side of (2) 
is the energy tensor for the photon fluid driving the rocket. To drive the rocket 
by gravitational waves we must get rid of this energy tensor, so that we solve 
the vacuum equations 

R''' = (3) 

Accordingly, we try to add terms to (1) so that the right-hand side of (2) dis- 
appears. 

The quantities on the right-hand side of (2) are 

r = SI, ^i = r-^{2M -6Mf cose); (4) 

where an overdot means differentiation with respect to u. In (4) the expressions 
2M and ~QMf cos 9 refer respectively to the loss of mass and the change of 
momentum by the rocket because of the emission of photons. We wish to add 
terms to the metric (1) which cancel these, so eliminating the photons which 
power the photon rocket. From the theory of gravitational radiation we have 
clues as to how to do this. We know from many investigations that a source loses 
mass by the emission of quadrupolc waves [14], and that it gains momentum 
from recoil when it emits quadrupolc and octupole waves [15]-[19]. This suggests 
that the terms we must add to (1) arc those describing quadrupolc and octupole 
gravitational waves. If we do this correctly we should get equations of motion 
for the rocket in terms of the gravitational radiation it emits: M and / will no 
longer be arbitrary functions. 

We shall add the quadrupolc and octupole wave terms in the linear 
approximation, but the cancellation of T**^ in (2) will occur in the quadrupole- 
quadrupole and quadrupole-octupole approximations. 

In Section 2 we describe the metric and approximation method which 
we use, and in Section 3 we explain the method of solving the field equations. 
Sections 4-7 outline the solutions of various approximation steps, the main re- 
sults appearing in equations (31) and (34). There is a brief concluding Section 



2 The approximation method 

We assume axial symmetry and take the m.etric in the form 

ds^ = -r^{Bde^ + Csin^ ed<j)^) + Ddu^ + 2Fdrdu + 2rGdedu; C = B"\ (5) 

where B, C, D, F, G are functions of r, 9, u. This metric, which is a form of that 
due to Bondi [20], is particularly suitable for describing outgoing radiation. It 
will be noticed that (1) is of the form (5). 

We use the double-series approximation method which has been used 
previously in conjunction with the metric (5) to study the emission of gravita- 
tional waves from isolated sources [19], so we assume that the solution we are 
looking for can be expanded in powers of constant parameters to, the initial 
mass of the rocket, and a, some convenient length associated with it. 

We envisage some internal vibrations of the rocket which generate vary- 
ing quadrupole and octupole moments Q{u) and 0{u) about the symmetry axis, 
given by 

Q = ma h(u), O = ma k{u) (6) 

where h and k are smooth functions which vanish except for a finite interval 
< u < Ml. There may, of course, be multipolc moments higher than the 
octupole, but we shall not be concerned with them here. 

The expansion of the metric functions takes the form 

oo oc 

ffjfe = X! X! "^^"^ 9tk , (7) 

p=0 s=0 

and we assume that M and / in (1) can also be expanded in powers of tti and 
a: 

p— 1 a— p— 1 a— 

(ps) (ps) (ps) (ps) (ps) 

In the above the Qik, M , f are independent of to. and a, and M , / are 
functions of u only. In fact, since the starting value of M is to. and decrements 

to the mass appear at higher approximations, we have M — I- 

Introducing the notation of (5) we write the metric coefficients as 

oo oo / ^ 

-.922 = r2B = r2[l + ;^^TOPa^ B], (9) 

oc oc , -. 

-333 = r^ sin^ 0C = r^ sin^ 6'[1 + ^ ^ mPfl" C], (10) 

p=l s=0 

oc oc , X 

544 = £> = 1 + E E ""^"^ ^ ' (11) 



OO CO / -. 

314 = F=l + ^5]mPa^ F, (12) 

OO OO , , 

324 = rG = r^^mPa^ G • (13) 

p=l s=0 

These summations ensure that if the mass parameter m is put zero we have flat 
spacetime in the form 

ds^ = ~r^{de^ + sin^ 9d(j)^) + 2drdu + du^. (14) 

We emphasise that we are looking for a metric of form (1) which satisfies 
the vacuum equations (3). There is now no energy tensor such as occurred in 
the photon rocket. 

If we substitute (7) into the field equations (3) we can separate out the 
coefficients of mPa^ . Equating these to zero, we obtain a doubly infinite set of 
second order differential equations. The set which is the coefficient of mPa^ will 
be called the (ps) approximation. The (00) approximation vanishes. The (Is) 

approximations arc linear and homogeneous in the g , and it is here that we 
shall insert the quadrupole and octupolc wave terms which are of order (12) 
and (13) respectively. The non- linear terms are those with p > 1. For further 
description of the double-series method see [19]. 

For the reason described in Section 1, we expect changes of mass to 
appear at the quadrupole-quadrupole, i.e. the (12)x(12) or (24) approxima- 
tion. Changes in momentum should come from the quadrupole-octupole, i.e. 
(12)x(13) or (25) approximation; the acceleration engendered by this change 
will be of order (15). Wc accordingly take the initial terms in (8) to be 

(24) 

M = m + rrCa^ M , (15) 

. (15) 
/ = ma^ f . (16) 



3 Solution of the field equations 

The (ps) approximation described above consists of seven equations of the form 

( 7JS ) ( QT' J 

^lm{gik)^'^lm{grk), (17) 

where the left hand side is linear in the gik (and their derivatives) , and the right 

(gr) 

hand side is non-linear in the gik , [q < p,r < s) and is known from previous 
approximation steps. For instance, in the (24) approximation, which we shall 



consider in Section 6, the terms on the right hand side could come from the 
following combinations 

(10) (14) (11) (13) (12) (12) 

9ik X gik , Qik X gik , gtk x gik (1»J 

(and their derivatives). 

In Appendix I of [19] were given the explicit forms of the left hand sides 
of (17), together with an algorithm for its solution. Here we give only a summary 
of the results. Where there is no ambiguity we shall drop the superscripts (ps) 
to save printing. 

The solution of the (ps) approximation depends on an inhomogcneous 
pseudo wave equation for D: 

U'D = Dii - 2Di4 + 2r-^(Di + D4) + r-^{D22 + -D2 cot e) = n + 2r"\4(u, 6) 

(19) 
where subscripts 1,2,4 mean differentiation with respect to r, 9, u, where 11 con- 
sists of non-linear terms of order (ps) known from previous approximations and 

(ps) 

also terms in F which is known by a simple integration of one of the [ps) 

(ps) 

equations. x(u, 0) is an arbitrary function of integration, and such a function 
arises at each approximation step; it enters not only the equation for D, but 
also, directly or indirectly, the equations for B and G. It plays an essential part 
in the construction of physically sensible solutions. 

(ps) (ps) (ps) (ps) 

The algorithm in [19] shows how to calculate B , C , G once D has 
been got from (19). In obtaining them one must ensure that there is no physical 
singularity on the rotation axis 9 = 0,0 = tt. Sufficient conditions for this are 
that 

B/si-a^e, C/sin^e, D, F, G/sinO (20) 

shall be of class G^ near sin 9 = [20] . 

4 The (Is) approximations 

The (10) approximation is simply the Schwarzschild solution in Bondi coordi- 
nates. In the notation we have described it is (omitting the (10) superscript) 

D = -2r-^, B = G = F = G = 0. (21) 

It was shown in [19] that the (11) approximation allows no non-singular 
solution except motion with constant velocity, which we shall put zero, so abol- 
ishing the (11) terms. This corresponds with the known result that there are 
known linear dipole gravitational waves. However, a non-linear dipole wave does 
appear at the (25) approximation, as we shall see in Section 7. 



The (12) and (13) approximations, corresponding to quadrupole and 
octupole waves were given in [19]: 

(12) (12) 1 

B = -C = -siT?e{r-^h + r-^h), (22) 

(12) 

D ^ -P2i2r-^h + 2r-^h + r-^h), (23) 

(12) 

F = 0, (24) 

^G = -P^{2r-^h-4:r-^h-3r'^h). (25) 



(13) (13) 1 , 1 ■■■ ,• 4 

B^~C = — sin2 6lcos6'(2r"n +10r"^fc + 15r"^fc, (26) 

D = -\p3{2r-^k+Ar-^k + 5r-^ + 3r-'^k), (27) 

(13) ^ ^ 

F = 0, (28) 

(13) 1 

G = —P^{2r'^ k -8r-^k - ISr^^fc _ Ur-^k). (29) 

oD 

Here P2(cos6'), P3(cos6') are Legendre polynomials, / means d/dd and h,k are 
the functions introduced in (6). 

The only other (Is), i.e. linear, approximation which enters the work 

(15) 

is the (15) one, which comes in because of / (u). We see from (1) and (16) 
that the (15) terms in the metric are (dropping the superscripts) 

B = C = F = 0, D = -2frcose, G ^ -frsmO. (30) 

If we substitute these expressions into (17) (with right-hand side zero because 
we are considering a linear approximation) we find it satisfied identically.Fl 

5 The second approximations 

Our methodology in this work, like that in [19] and [21], is to introduce the 
multipole waves in the linear approximation only. That means we do not add 
any further solutions of the homogeneous equations (i.e. (19) with vp/m = 0) in 
the higher approximations. This paradigm is ambiguous, because the solution 
of the non-linear equations (i.e. (19) with ^/i^ ^ 0) is not unique. This non- 
uniqueness of the double series method has been discussed in [21]. It applies 

^The terms (30) will give rise to a (2 10) approximation, as expected because of the ap- 
pearance of /^ in (1). This too vanishes identically. The satisfaction of the (15) and (2 10) 
approximations occurs because (1) with M = is flat. 



to other approximation methods which are governed by a set of equations such 
as (19), since at any approximation stage one can always add a solution of the 
homogeneous equations. Fortunately the non- uniqueness does not trouble us in 
the work described in this paper. 

Turning to the individual (2s) steps, we find there are no non-linear 

terms of order (20) or (21), so we choose gik and gik both zero. There are 
terms of order (22) and (23) but, as shown in [19], they refer to wave tails, 
which die away gradually to zero, and have no effect on the motion of the 
rocket. 

The first step of interest here is the (24) one, to which we now turn. 

6 The (24) approximation 

The possible non-linear contributions in the general (24) approximation are 

those in equation (18). However, in our work we are taking gik equal to zero, 
and inserting no (14) multipole, so the only non-linear terms we need are those 

(12) (12) 

corresponding to (7ifc x gik- This case was completely solved in [21]; the solution 
is lengthy and will not be reproduced here. It is free from singularities on the 
rotation axis, and consists of oscillating terms which cease when the internal 
motions cease, together with certain terms showing secular change. The latter 
can be removed by coordinate transformations, as was shown in [19], except 

(24) 

for one term in D which refers to a change in the Schwarzschild mass of the 
rocket, corresponding to the transmission of energy by the gravitational waves. 

(24) 

This is to be identified with M in equation (15): 

(24) 1 /•« ...2 

M-~yJ_^hdu. (31) 

Eqn (31) gives the mass-loss of the rocket in terms of its quadrupole oscillations. 
It agrees with emission of energy according to the quadrupole formula. 



7 The (25) approximation 



In solving the (25) approximation two conditions will be imposed. First, there 
must be no singularities on the rotation axis. Sufficient conditions for this were 
stated in equation (20). In the solution we are considering these can be satisfied 
by suitable choice of arbitrary functions, such as xi''^^ ^)i arising in the algorithm 
for solving each (ps) approximation. In general the question of whether x can 
be chosen to avoid singularities along the axis of symmetry depends on the 
satisfaction of a particular relationship between the non- linear terms at each 



order of the approximation. It has been shown that this relationship can be 
satisfied up to the fourth order in m and a [22]. 

Secondly, we must ensure that the rocket remains at the origin r = 
so that we know that its acceleration is /. This will be so if its dipole moment 
about the origin vanishes at all time. We shall assume that the dipole moment 
enters the (25) approximation as the coefficient of r~^ cos 6 in the expansion of 

(25) 

£) , as it would into the expansion of the Newtonian gravitational potential. 
It will be noted that the generating equation (19), in the linear approximation 
and the static case, (in which 11 and X4 are zero) reduces to Laplace's equation 
with axial symmetry, and has a solution fir^^ cos 6 for a dipole /i. 

The non-linear terms in the (25) approximation come from the following 
combinations 

(10) (15) (12) (13) 

gik X gik, gik X gik . 
The former of these consists of a single expression 

(15) 

—Grri^a^ f cos9/r^ 

arising in the Ricci tensor component R^^; in fact it corresponds to part of (4) 
which we wish to cancel, as explained in Section 1. The second combination 
arises from the quadrupole-octupole interaction. 

The (25) approximation has been completely solved by one of us (MSP). 

(25) 

Unfortunately the solution is extremely long (the expression for D alone con- 
tains 66 terms) so we do not print it here, but we will gladly send a copy to 
anybody who asks for it. The following is a brief summary. 

(25) 

D can be expressed exactly in a finite series 

(25) ^ 

i? = ^r-" A(u,0). 

n=l 

(25) 

The arbitrary functions which arise in solving the (25) step, including X (u, 0) 

(25) 

which enters the pseudo wave equation (19), must be chosen so that G and 

(25) 

B are regular on the axis. When this is done, we find 



(25) , 2 r ^^'"'^ 

D = cose[r-'{-—Y - 6 / (fix) dx) 

105 Jn 



j-u i-x (15) 2 2 1 

+r-\-2 I I ifiy) dydx + —Z + —E - —p{u)q{u)) + A{r, u)] 



JO 



+P3X{r,u) + P5W{r,u). (32) 



In this 

p{x) := (fh{x)/dx'^, q{x) ■- <fk{x)/dx^ , 

Y= r d_P^d_q^dx, Z^ r Y{x)dx, E^ T ^,(x)dx, (33) 



J r\ ax ax J r\ J r\ ax 

A denotes regular transient terms of order r~^ and higher, and W, X are finite 
series in r~^ with regular transient coefficients which need not concern us here. 
We now impose the second condition referred to above, namely that the 
coefficient of r^^ cos9 shall vanish. Differentiating this twice with respect to u 
we have 

(15) 1 

630 
where p, q are functions of u. So using equation (16) and introducing units of 
customary dimensions 

f - ggo^ (2m - ipq + 3qp), (34) 

G being here the gravitational constant. This gives the acceleration of the rocket 
in terms of the rates of change of the quadrupole and octupole moments induced 
by the rocket's internal motions correct to order (15) in the expansion parame- 
ters m and a. 

(25) ^ 

in D is 

cos 9 

-:r^(pq-qp), 
35r 

which represents a dipole gravitational wave. Thus although a dipole wave is 
forbidden in the linear approximation, one arises at the (25) stage. 

(25) (25) 

The solutions for B and G obtained by the algorithm in [19] con- 
tain secular terms involving integrals like those occurring in equation (32). 
However, they can be removed by a coordinate transformation, namely the 
a-transformation given in [20], which was also used in [19], eqn (10.14). 

The acceleration lasts while there are internal motions, i.e. from u = 
till u — ui. Integrating the acceleration (34) between these limits we have the 
velocity V acquired, assuming the rocket was at rest at m = 0: 

^, mGa"^ n , 
^=315^70 """ ' 
agreeing with a result obtained in [19]. 



8 Conclusion 

The main results of the paper are equations (31) and (34), giving respectively 
the mass-loss and acceleration of the rocket in terms of its internal motions 
represented in our work by variations of the quadrupolc and octupole moments. 
We have considered the following approximations 

(00), (10), (12), (13), (15), (24), (25), (2 10) 

and referred briefly to the wave-tail approximations (22) and (23) which do not 
affect the equations of mass- loss and motion of the rocket. There will of course 
be other approximation steps, even if we confine ourselves to the quadrupole 
and octupole oscillations. For example, the (26) approximation arising from 
(13) X (13) terms will contribute to the mass-loss. 
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